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Using numerical simulations, we demonstrate that superconductors containing periodic pinning 
arrays which have been diluted through random removal of a fraction of the pins have pronounced 
commensurability effects at the same magnetic field strength as undiluted pinning arrays. The 
commensuration can occur at fields significantly higher than the matching field, produces much 
greater critical current enhancement than a random pinning arrangement due to suppression of 
vortex channeling, and persists for arrays with up to 90% dilution. These results suggest that 
diluted periodic pinning arrays may be a promising geometry to increase the critical current in 
superconductors over a wide magnetic field range. 

PACS numbers: 74.25.Qt 



I. INTRODUCTION 

One of the most important issues for applications of 
type-II superconductors is achieving the highest possi- 
ble critical current. This requires preventing the depin- 
ning and motion of the superconducting vortices that are 
present in the sample. A promising approach is the use 
of artificial pinning arrays, which have been the focus of 
extensive recent studies. Pinning arrays with triangu- 
lar, square, and rectangular geometries have been fabri- 
cated in superconductors using either microholes or blind 
holesii^ii^i^i^i^ or arrays of magnetic dots^i^. The resulting 
critical currents are significantly enhanced at the match- 
ing magnetic field where the number of vortices equals 
the number of pinning sites, as well as at higher fields 
nB^, with n an integer. At the matching fields, peaks 
or anomalies in the critical current occur and the vor- 
tices form highly ordered commensurate lattices which 
are free of topological defects^i^i^i^i^ii^iiiii^ii^. Although 
it might appear that the best method for increasing the 
overall critical current would be to increase the pinning 
density and thus increase B^, this technique fails since 
adding too many pins to the sample can degrade the over- 
all superconducting properties by lowering Tc and Hc2- 
Thus, an open question is how the critical current can be 
maximized using the smallest possible number of pinning 
sites. Since completely periodic pinning arrays have been 
shown to enhance pinning, one can ask whether other pin- 
ning arrangements such as quasiordered or semiordered 
arrays could be even more effective at pinning vortices. In 
recent simulationsi^ and experiment a^^d^ , it was demon- 
strated that a quasicrystalline pinning array produced an 
enhancement of the critical current compared to a ran- 
dom pinning array for fields below and up to B^. The 
enhancement disappears for fields above B^, however. 

In this work we demonstrate a new type of commensu- 
rability effect that occurs in periodic pinning arrays that 
have been diluted by randomly removing a fraction of the 
pinning sites. In an undiluted array, the commensurabil- 
ity effects occur at fields ni?^. In a diluted array, i?^ is 
reduced; however, the pinning array still retains corre- 



lations which are associated with the periodicity of the 
original undiluted pinning array. In this case, in addition 
to commensuration effects at B^, we observe noticeable 
commensuration effects at the higher field B^, the match- 
ing field for the undiluted pinning array, even though 
many of the pins that would have been present in such 
an array are missing. Strong peaks in the critical cur- 
rent associated with well-ordered vortex configurations 
appear at nB^. This commensurability effect is remark- 
ably robust and still appears in arrays which have up to 
90% of the pinning sites removed, so that — lOB^. In 
samples with equal numbers of pinning sites, diluted pe- 
riodic pinning arrays produce considerable enhancement 
of the critical current for fields well above B^ compared 
to random pinning arrays. The diluted periodic pinning 
arrays also have a reduced amount of one-dimensional 
vortex channeling above B^ compared to undiluted peri- 
odic pinning arrays. Such easy vortex channeling is what 
leads to the reduction of the critical current above B^ in 
the undiluted arrays. Our results should also apply to 
related systems such as colloids interacting with periodic 
substratesi^^ and vortices in Bose-Einstein condensates 
interacting with optical trapsi^. 



II. SIMULATION 

We simulate a two-dimensional system of size x Ly 
with periodic boundary conditions in the x and y direc- 
tions containing iV„ vortices and Np pinning sites. The 
vortices are modeled as point particles where the dynam- 
ics of vortex i is governed by the following equation of 
motion: 

V^^Fr+F^ + F^. (1) 

The damping constant rj — 0§d/27r^^pjv, where d is the 
sample thickness, ^ is the coherence length, is the 
normal-state resistivity, and (pQ = h/2e is the elementary 
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flux quantum. The vortex-vortex interaction force is 



F™ -^A„/oifi(i?,,/A)R,, 



(2) 



where Ki is the modified Bessel function, Ay is an inter- 
action prefactor which is set to 1 unless otherwise noted, 
A is the London penetration depth, /o = (jy^ / 2i: ^qX^ , 



is the position of vortex Rij 



IR," — Ri 1, and 



(Ri — Ilj)/Rij. Forces are measured in units of 



R 
R 

/o and distances in units of A. The vortex density B — 
Ny/LxLy. The vortex-vortex interaction force falls off 
sufficiently rapidly that a cutoff at Rij = 6A is imposed 
for computational efficiency and a short range cutoff at 
Rij = O.IA is also imposed to avoid a force divergence. 
The Np pinning sites are modeled as attractive parabolic 
potential traps of radius rp = 0.2A and strength /p = 

0.25/0, with F^^" = (/p/rp)i?.fce((rp - i?,fc)/A)f|^^ Here 

R^^^ is the location of pinning site fc, Riu — |Ri — Rif^l, 

Rifc — (Ri — R[,^'')/i?ifc, and O is the Heaviside step func- 
tion. The pinning density is ~ Np/{LxLy). The pins 
are placed in a triangular array with matching field B^ 
that has had a fraction Pd of the pins removed, so that 
B^ = {1- Pd)B*^. The Lorentz driving force F^ = F^x, 
assumed uniform on all vortices, is generated by and per- 
pendicular to an externally applied current J. The ini- 
tial vortex positions are obtained by simulated annealing 
from a high temperature. After annealing, the tempera- 
ture is set to zero and the external drive is applied in 
increments of 2 x 10~^/o every 5 x 10^ simulation time 
steps. We measure the average vortex velocity in the di- 
rection of the drive, ly) = '^i'^i where is the 
velocity of vortex i, and define the critical depinning force 
fc as the value of F'" at which {V) = 0.01. We have found 
that for slower sweep rates of the driving force, there is 
no change in the measured depinning force. Finally, we 
note that our model should be valid for stiff vortices in 
three-dimensional superconductors interacting with ar- 
rays of columnar defects. For a strictly two-dimensional 
superconductor with periodic pinning arrays the vortex- 
vortex interaction is modified from a Bessel function to 
ln{r); however, our previous simulations have indicated 
that either interaction produces very similar commensu- 
rability eSectsS^^^^i^^. 
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FIG. 1: Circles: Pinning site locations for two 18A x 18A 
samples with the same number of pinning sites Np = 168. 
(a) A triangular pinning array with no dilution, Pd — 0, with 
B4, = 0.52/A^ and = 0.52/A^. (b) A triangular pinning 
array that has been diluted at Pd — 0.5, with — 0.52/A^ 
and = 1.04/A^ 
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FIG. 2: The critical depinning force fc/fo versus B/B^ for 
systems with B^ — 0.52/A^. The actual matching field B^ 
and the matching field for an equivalent undiluted array 
are labeled, (a) /c//o for an undiluted pinning array with 
Pd = G (heavy line), a pinning array with Pd = 0.2 (light 
line), and a random pinning arrangement (dashed line), (b) 
Pd = (heavy line) and Pd = 0.4 (light line), (c) Pd = 
(heavy line) and Pd = 0.5 (light line), (d) Pd = (heavy line) 
and Pd = 0.65 (light line). 



III. INCREASING DILUTION 

A. Fixed Pinning Density 

In Fig. 1(a) we illustrate the positions of the pinning 
sites for an 18A x 18A system with no dilution, Pd = 0. 
The pinning sites are placed in a triangular array at a 
pinning density of = 0.52/A^, giving B^ = B'^ = 
0.52/A^. Figure 1(b) shows a system with B^ — 0.52/A^ 
and dilution Pd = 0.5, so that B*. = 1.04/A2. Here, 



half of the pinning sites that would have formed a denser 
triangular array have been removed randomly to give the 
same matching field as in Fig. 1(a). 

To illustrate the effect of the dilution, in Fig. 2(a) we 
plot the critical depinning force fc/fo vs B/ B^ for three 
systems with B^j, = 0.52/A^. Two of the samples contain 
triangular pinning arrays at different dilutions, Pd = 
[as in Fig. 1(a)] and Pd = 0.2. The third sample contains 
randomly placed pinning. For the undiluted triangular 
pinning array with Pd = 0, peaks in the depinning force 
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occur at integer values of B/B^. A submatching peak 
at B/B^ = 1/3 also appears in agreement with earlier 
studiesi^. The random pinning arrangement shows no 
commensurability peaks and has a lower critical depin- 
ning force than the periodic pinning array for all but the 
very lowest values oi B / B(j,. The diluted periodic pinning 
array has a pronounced peak in fc/fo at B / B^ = 1.25. 
This peak corresponds to the matching field i?J of the 
equivalent undiluted pinning array. A second peak in 
fc/fo occurs at B — 2B'^. There is a still a small peak 
in fc/fo s± B / B^ = 1 for the diluted array; however, we 
do not find a peak at B/B^ = 2. The diluted pinning 
array has a higher fc/fo than both the random pinning 
arrangement and the undiluted periodic pinning array for 
B/B^ > 1 over the range of B/B^ investigated here. In 
most of this range, fc/fo for the diluted periodic array is 
twice that of the undiluted periodic array. 

The peak in fc/fo at B = B'^ persists as the dilution 
increases. Figure 2(b) shows that a diluted pinning array 
with Pd = 0.4 has a small peak at B/B^ — 1 and a more 
prominent peak at B = B^. The critical dcpinning force 
is again higher in the diluted array than in the undiluted 
array for B / B^ ^ l-l- For a- diluted pinning array with 
Pd = 0.5, Fig. 2(c) indicates that the commensurability 

peak in fc/fo has shifted to B = B* = 2Bs. Here, the P^^' Positions of pinning sites (open circles) and vortices 



ee®®.®®.®® 
® ■ . ® ■ ® . . ■ ® 



■®.®®®®.®®®® 



(C) 



(d) 



"0 

pinning has become so dilute that there is no longer a 
noticeable enhancement of fc/fo at B/B^ = 1. When 
Pd = 0.65, as in Fig. 2(d), the peak in fc/fo shifts up to 
B ^ B'^ = 2.86i?0. As increases with increasing dilu- 
tion, the critical depinning force at lower fields B < B^ 
decreases. These results suggest that diluted periodic ar- 
rays may be useful for increasing the overall pinning force 
at higher fields, and that a peak in the critical current 
at a specific field Bp can be achieved using a significantly 
smaller number of pins than would be required to create 
an undiluted periodic pinning array with B^ — Bp. 

To understand the origin of the peak in fc/fo at _B^, 
we analyze the vortex configurations at B*^ for four dif- 
ferent values of Pd in Fig. 3. For an undiluted triangular 
pinning array at Pd — 0, as shown in Fig. 3(a), B*^ — B^, 
each pin captures exactly one vortex, and there are no in- 
terstitial vortices trapped in the regions between pinning 
sites. In Fig. 3(b), at B = B*^ — 1.25B^ in a system with 
Pd = 0.2, all of the pinning sites are occupied and the ex- 
cess interstitial vortices sit at the locations of the missing 
pinning sites, creating a triangular vortex lattice. As Pd 
increases, the increasing number of interstitial vortices 
present at -B = B'^ continue to occupy the locations of 
the missing pinning sites. This is illustrated for Pd = 0.5 
in Fig. 3(c) and Pd = 0.65 in Fig. 3(d). The high symme- 
try of the triangular vortex lattice configuration at the 
matching field B = B'^ causes the vortex-vortex inter- 
actions to cancel, and the depinning force is determined 
primarily by the maximum force of the pinning sites, fp. 
Defects in the vortex lattice appear just above and be- 
low B = creating asymmetrical vortex- vortex inter- 
actions. A vortex associated with a defect experiences 
an extra force contributing to depinning on the order of 



(black dots) for samples with 
2B^. (d) Pa 



0.52/A^ at B/S; = 1. (a) 



Pd = 0, S; = B^. (b) Pa = 0.2, B; = 1.255^. (c) Pd = 0.5, 



b; 



0.65, b; 
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F'"'"{a/\), where a is the vortex lattice constant, and fc is 
reduced. A similar disordering process occurs just above 
and below each of the higher matching fields^, although 
in this case the depinning force is not merely determined 
by the pinning force since a portion of the vortices are 
pinned indirectly in interstitial positions and experience 
a weaker effective pinning force. 

Previous work on wire network arraysiS showed that 
even if the network is partially disordered, spatial corre- 
lations remain present as indicated by the appearance of 
peaks in fc-space and can produce matching effects. In 
the diluted pinning arrays we consider here, even though 
a portion of the sites are removed there are still peaks in 
the reciprocal space corresponding to the original undi- 
luted array, so matching effects occur when the vortex 
lattice k spacing is the same as these pinning array k 
space peaks. We note that if an equivalent number of 
randomly placed pinning sites is used there is no match- 
ing effect since there is a ring rather than peaks in k 
space. 

In the case of the diluted periodic pinning array, the 
vortex configuration contains numerous topological de- 
fects at B — B^ for higher values of Pd when the pinning 
is strong enough that most of the vortices are trapped 
at the pinning sites. In contrast, at i? = B*^ the vor- 
tex lattice is free of defects, as illustrated in Fig. 3. For 
fields just above or below B = B^, vacancies and intersti- 
tials appear in the vortex lattice and lower the depinning 
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force. At higher fields B = nB'^, where n > 1 is an inte- 
ger, the vortex lattice is again ordered and a peak in the 
critical depinning force occurs, as shown in Fig. 2(a) at 

B = 2b;. 

The enhancement of fc/fo in the diluted periodic pin- 
ning arrays for B / B^j, > 1 at the nonmatching fields oc- 
curs due to a different mechanism. We first note that at 
B/Bcj, < 1, fc/fo is lower in the diluted pinning arrays 
than in the undiluted pinning arrays since it is possi- 
ble for some of the vortices to sit in interstitial sites in- 
stead of in pinning sites, even though not all of the pins 
have been filled. Interstitial vortices experience an effec- 
tive pinning force due to the interactions with the sur- 
rounding vortices. Since the interstitial vortices are more 
weakly pinned than vortices in pinning sites, the critical 
depinning force is lower than it would be for an undiluted 
array at B / B^ < 1, where every vortex sits in a pinning 
site. For B/B^ > 1, the initial motion of the vortices 
at depinning in an undiluted array occurs via channel- 
ing of the interstitial vortices between the pinning sites 
along an effective modulated one-dimensional potential 
created by the vortices located at the pinning sitesii. In 
a diluted periodic pinning array at i?^ < B < iJJ, any 
channels of interstitial vortices are interrupted by one or 
more pinning sites which serve to increase the depinning 
threshold of the entire channel by "jamming" free motion 
along the channel. Here, the one-dimensional potentials 
that allow for relatively free interstitial vortex motion 
do not form until B > B^ when the vortices begin to 
occupy positions that would be interstitial sites in the 
equivalent undiluted pinning array. This picture is con- 
firmed by our analysis of the vortex trajectories. The 
jamming effect eventually breaks down at high dilutions 
Pd when large regions devoid of pinning sites span the 
entire system and provide macroscopic channels for free 
vortex fiow. In Fig. 2(d) at Pa = 0.65, fc/fo at B/B^ > 1 
is lower than the fc/fo values obtained at lower P^ and 
similar fields. For higher values of Pd this effect becomes 
more pronounced until fc/fo eventually drops below the 
critical depinning current for the undiluted array. 



B. Decreasing Pinning Density 

In Fig. 4 we show the effect of gradually diluting a peri- 
odic pinning array with — 0.83/A^. The main peak in 
the critical depinning current always falls at B = B^, but 
as the dilution increases from = to = 0.9, the 
value of fc/fo decreases significantly. The pronounced 
peak in fc/fo persists even up to 90% dilution. There is 
a small peak at B = B^ for Pd = 0.25, but for higher dilu- 
tions we find no noticeable anomalies in fc/fo at P = B^. 
This result indicates that the commensurability effect at 
B = B1 is remarkably robust. 
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FIG. 4: fc/fo vs B/B*^, where B; = 0.83/A^ Top curve: 
an undiluted periodic pinning array with B^ — 0.83/A^. Re- 
maining curves, from top to bottom: Successive dilutions of 
the same array to Pd = 0.25, 0.5, 0.75, and 0.9. 
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FIG. 5: fc/ fp vs Av, the vortex- vortex interaction prefactor, 
for four samples with B = 0.52/A^. Circles: undiluted tri- 
angular array with B/B^ — 1 and B^ = 0.52/A^. Squares: 
diluted pinning array with Pd — 0.5 at B/B'^ — 1 where 
B^ = 0.52/A^ Triangles: undiluted array with B/B^ = 2.0 
and B^ — 0.26/A^. Diamonds: random array with B/B^ = 
2.0 and B^ = 0.26/A^ 

IV. EFFECT OF VORTEX LATTICE STIFFNESS 

We next consider the effect of the vortex lattice stiff- 
ness. This parameter can be changed artificially by ad- 
justing the value of the factor Ay in the vortex-vortex 
interaction term of Eq. (2). For dense random pin- 
ning where there are more pinning sites than vortices, 
Np > Ny , decreasing the strength of the vortex- vortex in- 
teractions by lowering Ay makes the vortex lattice softer 
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and the depinning force increases. In the limit oi Ay = 
the vortices respond as independent particles and the de- 
pinning force will equal fjr^'^- In the case where there 
are more vortices than pinning sites Ny > Np, the situ- 
ation can be more complicated since there are two types 
of vortices. One species is located at the pinning sites 
and the other species is in the interstitial sites. The in- 
terstitial vortices are not directly pinned by the pinning 
sites but are restrained through the vortex-vortex inter- 
action with pinned vortices. In this case it could be ex- 
pected that increasing the vortex-vortex interactions by 
raising would increase the depinning force. On the 
other hand, if the vortex- vortex interaction force is strong 
enough, the interstitial vortices can more readily push the 
pinned vortices out of the individual pinning sites. Stud- 
ies of very diluted random pinning arrays found that as 
the vortex-vortex interaction is increased, the depinning 
force initially increases linearly; however, for high enough 
vortex-vortex interactions the depinning curve begins to 
decrease again^^. The peak in the depinning force co- 
incides with a crossover from plastic depinning at low 
vortex-vortex interaction strengths to elastic depinning 
at high vortex-vortex interaction strengths. A similar 
behavior was observed for vortices in periodic pinning 
where there were more vortices than pinning sites^i. 

For the results presented in Fig. 2 and Fig. 4, the 
vortex- vortex interaction prefactor was set to A^ = 1.0 
and the depinning was plastic with the interstitial vor- 
tices depinning first. In Fig. 5 we plot fc/fp vs Ay at 
fixed B = 0.52/A^ for four different pinning geometries. 
In an undiluted pinning array with B/B^ = 1, the vor- 
tex lattice matches exactly with the pinning array and 
the vortex-vortex interactions completely cancel. As a 
result, the system responds like a single particle and 
fc = fp for all Ay, as shown by the circles in Fig. 5. 
When the same pinning array is diluted to = 0.5 and 
held at B/B'^ = 1, the squares in Fig. 5 illustrate that 
for 0.25 < Ay < 5.0 the depinning is plastic and fc/fp 
linearly increases with Ay. This behavior is agreement 
with earlier studies^i^^. For Ay > 5.0 the depinning is 
elastic and the entire lattice moves as a unit. In this elas- 
tic regime fc/fp first levels off and then decreases with 
increasing Ay as the force exerted on the pinned vortices 
by the interstitial vortices increases. This result indicates 
that fc/fp is maximized at the transition between plastic 
and elastic depinning. The filled triangles in Fig. 5 show 
the behavior of an undiluted array at B/B^ = 2 where 
the vortex lattice has a honeycomb structure rather than 
triangular symmetry^. Here the fc/fp versus Ay curve 
has a very similar trend to the Pd = 0.5 case, with a peak 
in fc/fp at the crossover between plastic and elastic de- 
pinning followed by a decrease in the depinning force for 
Ay > 5.0; however, fc/fp drops off more quickly for the 
undiluted pinning array at B/S^ = 2 than for the di- 
luted pinning array with Pd = 0.5. The vortex lattice is 
always triangular for the Pd = 0.5 diluted pinning array 
at B/B'^ — 1, so as Ai, increases there is no change in the 
vortex lattice structure and every pinning site remains 



occupied. For the undiluted pinning array at B/B^ = 2 
the vortex lattice experiences an increasing amount of 
distortion for Ay > 5.0 as the structure changes from a 
honeycomb arrangement in which all of the pinning sites 
are occupied to a triangular arrangement in which some 
pins are empty. This implies that for stiff vortex lattices, 
a relatively larger depinning force can be obtained from a 
diluted pinning array than from an undiluted pinning ar- 
ray with an equal number of pinning sites. The triangles 
in Fig. 5 show the behavior of a random pinning array 
at a pinning density of B^ — 0.26/A^ with B/B^ = 2. 
We find the same trend of a peak in fc/fp as a function 
of Ay] however, fc/fp for the random pinning array is 
always less than that of the periodic or diluted pinning 
arrays. We note that a peak in fc/fp still appears for 
vortex densities that are well away from commensura- 
tion due to the competition between the pinning energy 
and clastic vortex lattice energy. 

One aspect we did not explore in this work is the ef- 
fect of melting and thermal fluctuations on the diluted 
pinning arrays. Melting would be relevant in strongly 
layered superconductors with periodic arrays of colum- 
nar defects. If a diluted periodic pinning array could be 
fabricated in a layered superconductor, it is likely that 
a rich variety of thermally induced effects would appear. 
One system that should have similar behavior is a sample 
with a low density of randomly placed columnar defects 
in which the number of vortices is much greater than 
the number of pinning sites??. Theoretical work sug- 
gests that such a system can undergo a Bragg-Bose glass 
transitioii^ii^. In the case where the pinning is periodic, 
it would be interesting to determine whether the melting 
transition shifts upward or downward at B / B^ ~ 1 com- 
pared to away from B/B'^ = 1. Additionally, at inter- 
mediate dilution it may be possible to observe two-stage 
melting or a nanoliquid. Such a nanoliquid has been ob- 
served in experiments with patterned regions of columnar 
defects^^. 



V. SUMMARY 

In summary, we have demonstrated that commensu- 
ration effects can occur in superconductors with diluted 
periodic pinning arrays. Pronounced peaks in the critical 
depinning current occur at magnetic fields corresponding 
to the matching field for the undiluted array, B*^, which 
may be significantly higher than the matching field B^f, for 
the pinning sites that are actually present in the sample. 
The commensurability effect is associated with the for- 
mation of ordered vortex lattice arrangements containing 
no topological defects. This effect arises since the diluted 
pinning array has correlations in reciprocal space at the 
same values of k as the undiluted pinning array, enabling 
a matching effect to occur for the diluted array when the 
vortex density matches the density of the undiluted array. 
The effect is remarkably robust, and peaks in the critical 
depinning force persist in pinning arrays that have been 
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diluted up to 90%. For small dilutions, a peak can also 
appear at the true matching field B^. In samples with 
equal numbers of pinning sites, diluted periodic arrays 
have smaller critical depinning currents than undiluted 
arrays at fields below B^, but produce a significant en- 
hancement of the critical current at fields above B^ com- 
pared to both undiluted arrays and random pinning ar- 
rangements. This enhancement for B < B^ is due to the 
suppression of easy channels of one-dimensional motion 
that occur in the undiluted arrays. We have also exam- 
ined the effect of the vortex lattice stiffness on depinning. 
For B > Bcf,, the depinning is plastic in soft lattices and 
the interstitial vortices begin to move before the pinned 
vortices, while for stiff lattices the depinning is elastic 
and the entire lattice depins as a unit. In the plastic de- 
pinning regime which appears at small lattice stiffness, 



the depinning force increases with increasing lattice stiff- 
ness since the interstitial vortices become more strongly 
caged by the vortices at the pinning sites. In the elastic 
depinning regime the depinning force decreases with in- 
creasing lattice stiffness. A peak in the depinning force 
occurs at the crossover between plastic and elastic depin- 
ning. For stiffer lattices the diluted pinning arrays show 
more pronounced matching effects. Our results suggest 
that diluted periodic pinning arrays may be useful for 
enhancing the critical current at high fields in systems 
where only a limited number of pinning sites can be in- 
troduced. 
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